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Modified Newtonian dynamics (MOND) is an empirical theory originally proposed to explain the 
rotation curves of spiral galaxies by modifying the gravitational acceleration, rather than by invoking 
dark matter. Here,we set constraints on MOND using an up-to-date compilation of kinematic tracers 
of the Milky Way and a comprehensive collection of morphologies of the baryonic component in the 
Galaxy. In particular, we find that the so-called “standard” interpolating function cannot explain 
at the same time the rotation curve of the Milky Way and that of external galaxies for any of 
the baryonic models studied, while the so-called “simple” interpolating function can for a subset 
of models. Upcoming astronomical observations will refine our knowledge on the morphology of 
baryons and will ultimately confirm or rule out the validity of MOND in the Milky Way. We also 
present constraints on MOND-like theories without making any assumptions on the interpolating 
function. 


If Newtonian gravity holds, the rotation curve of the 
Milky Way cannot be explained by visible (baryonic) 
matter only, therefore providing evidence for dark mat¬ 
ter in the Galaxy, from its outskirts (e.g. mi) down 
to inside the solar circle [1]. In this paper we derive 
constraints on modifications of Newtonian gravity based 
on the same data set as used in Ref. [3|- We focus our 
attention on modified Newtonian dynamics (MOND), 
originally proposed by Milgrom m (see also (8HT21). 
MOND postulates that below a characteristic value do 
the acceleration a is modified with respect to that pre¬ 
dicted by Newton’s law aiv by the introduction of an 
interpolating function /i such that 13 El 

'‘W° = a "’ (1) 

where fi(x) ~ x for i« 1 and pi{x) ~ 1 for x 1 (for 
modern theories of MOND, see El)- The theory does not 
predict the value of do nor a specific functional form for 
/x from first principles, so different proposals have been 
made in the literature (see e.g. [13 EU). Two widely 
discussed functional forms for /x are the “standard” in¬ 
terpolating function, 

X 

Rstd(x) = , =, ( 2 ) 

Vl + ar 

and the “simple” interpolating function, 

X 

Tsim(x) = Z “ • (3) 

1 + X 

The limit of large accelerations, i.e. a > ao or x > 
1, is strongly constrained by solar system tests, where 


MOND is bound to recover Newtonian gravity. The limit 
of small accelerations, i.e. a do or x <C 1, is usually 
invoked in MOND studies to highlight how the observed 
flat rotation curves of spiral galaxies can be reproduced 
in the absence of dark matter. 

In this paper, we set out to test the most common 
MOND scenarios with the latest data on the rotation 
curve of the Galaxy. We use a comprehensive compilation 
of kinematic tracers of the Milky Way and a state-of-the- 
art modelling of the baryons, both presented in Ref. [2]. 
Our results and analysis technique are complementary 
to previous MOND analyses of Milky Way data [131116] 
(see also [MU for constraints from the vertical force in 
the discs of galaxies, including our own). The quantita¬ 
tive study of alternative formulations of MOND or other 
theories of modified gravity is left for future work. 

There are two observational inputs required to test 
MOND in our Galaxy: the observed gravitational accel¬ 
eration a and that predicted by Newtonian gravity a at. 
These are the two key ingredients of our analysis. The ac¬ 
celeration a is obtained from the rotation curve through 
a = Ru ; 2 , where R is the Galactocentric radius and the 
angular circular velocity co c is set by the compilation of 
rotation curve measurements of Ref. [4] . This is a com¬ 
prehensive compilation of kinematic tracers optimised to 
R = 3 — 20 kpc and that includes an unprecedented set 
of published data on gas kinematics, star kinematics and 
masers. A detailed description of the data used and their 
treatment can be found in the supplementary information 
of Ref. [2] (cf. Tab. SI therein). The Newtonian accel¬ 
eration ajv = Rw 2 is set by the three-dimensional den¬ 
sity distribution of baryons in our Galaxy for which we 
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FIG. 1: The rotation curve of our Galaxy. The red data points indicate the angular circular velocity of all tracers in the 
rotation curve compilation, while the grey band brackets the contribution of all baryonic models under the assumption of 
Newtonian gravity. For further details, see Ref. [4]. We have adopted here Ro = 8 kpc, Vo = 230km/s and (U, V, IT) 0 = 
111.10.12.24. 7.251 kin/s EDI. 


adopt the survey of observation-based models presented 
in Ref. [4]. This survey includes seven alternative mor¬ 
phologies for the stellar bulge [IHTThTi , five for the stellar 
disc m mm and two for the gas [3TJ 152] , while the 
normalisation is set by microlensing observations for the 
bulge [33] [34], the local total stellar surface density for 
the disc [THj and the CO-to-H 2 factor for the gas ?>5 361. 
For further details, please refer to the supplementary in¬ 
formation of Ref. [1]. We use the morphologies of bulge, 
disc and gas in all 70 possible combinations, thus includ¬ 
ing all configurations of the baryonic component of the 
Galaxy present in the literature. In this way, we bracket 
the current uncertainty due to baryonic modelling and 
our conclusions do not rely on a specific model of the 
visible Galaxy but are solid against baryonic systematics. 
We adopt a distance to the Galactic centre Rq = 8 kpc, 
a local circular velocity Vo = 230km/s and the pecu¬ 
liar solar motion (t/, V, W) 0 = (11.10,12.24, 7.25) krn/s 
[2U] . The impact of varying these Galactic parameters is 
quantified later on, showing that our conclusions are ro¬ 
bust against current uncertainties. Only kinematic trac¬ 
ers with Galactocentric radii above R cut = 2.5 kpc are 
considered, amounting to N = 2686, 2687, 2715 individ¬ 
ual measurements for R 0 = 7.98, 8, 8.68 kpc, respectively 
(cf. below for the adopted uncertainty on R 0 ). The ef¬ 
fect of pushing the radius cut to R cut = 4.5 kpc is also 
quantified in order to avoid any influence of the bar (see 
e.g. Ref. (37]); in this case, N = 2159, 2162, 2267 for 
R 0 = 7.98, 8, 8.68 kpc, respectively. 


Fig-0 shows w c , i.e. the observed rotation curve, and 
Wb, i.e. the rotation curve expected from baryons un¬ 
der Newtonian gravity. We are now in place to com¬ 
pare these two observation-based quantities and set con¬ 
straints on MOND by assessing the discrepancy from the 
Newtonian scenario. The rotation curve predicted by a 
given MOND theory with fixed fi(x) and a o is obtained 
by solving Eq. 0 for a = -Rw^ond- We then quantify 
the goodness-of-fit of the obtained w mon d(R) curve using 
a two-dimensional chi-square against the rotation curve 
data w c (R). The use of this test statistic is mainly mo¬ 
tivated by the sizeable uncertainties on both u) c and R 
(see Ref. [1] for a technical discussion). 

The reduced chi-square as a function of the assumed 
cio is shown in the upper panels of Fig. [2] for each bary¬ 
onic model for n = /z st d (left) and /r = /z s ; m (right). The 
results for what we shall call fiducial baryonic model I 
m m ini are highlighted in black. Also shown are 
the values of a o favoured by the observation of rotation 
curves in external galaxies, namely ao = (1.27 ± 0.30) x 
10~ lo m/s 2 for fi s td and ao = (1.22 ± 0.33) x 10” lo m/s 2 
for n s ; m [38], as well as the 5cr exclusion line to guide the 
eye. In the case of the standard interpolating function, 
an acceptable best fit is obtained for most baryonic con¬ 
figurations. However, even in those cases, the preferred 
values of ao lie in the range (2.5 — 5.0) x 10 _lo m/s 2 , 
clearly above the typical values found in external galax¬ 
ies, which disfavours the standard interpolating function 
as being able to accommodate the rotation curve of ex- 
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FIG. 2: Fitting MOND to the rotation curve of our Galaxy. The left (right) panels show the reduced chi-square of the MOND 
scenario with the standard (simple) interpolating function for different values of ao- The upper panels convey the results for 
all baryonic models adopting R 0 = 8kpc, Vo = 230km/s and (U, V, W) 0 = (11.10,12.24, 7.25) km/s [20]. The black line in each 
upper panel marks the reduced chi-square obtained for fiducial baryonic model I 1191 j2ll [3lj . The bottom panels correspond 
to fiducial baryonic model I [191 1211131] with different combinations of Galactic parameters and radius cut. The range of a o 
found for each interpolating function from the rotation curves of external galaxies [3S] is encompassed by the vertical dashed 
lines, while the thick red line indicates the reduced chi-square corresponding to a 5a exclusion. 


ternal galaxies and that of the Milky Way at the same 
time. Instead, in the case of the simple interpolating 
function, the best-fit values of ao drop significantly to 
(1.0 — 3.5) x 10 -10 m/s 2 , in line with the range inferred 
from external galaxies for a subset of baryonic models. 
Nevertheless, most baryonic models still prefer some¬ 
what large values of ao- These results are roughly in line 
with previous fits of MOND to Milky Way data [El El, 
but here we have expanded upon such analyses by using 
a comprehensive compilation of the available kinematic 
data and a wide range of baryonic models. 

So far we have kept fixed Ro = 8kpc, vq = 230 km/s 
and Vq = 12.24 km/s. These values lie well within 
the ranges encompassing most current measurements, 
namely Rq = 8.0 ± 0.5kpc vq = 230 ± 20km/s 

02H47] and F 0 = 5.25 - 26 km/s [23 S3 SZ1 HE] ■ It is 
important to understand how the uncertainties on l?o, vq 
and Vq affect the results. For concreteness, let us focus 
on two specific measurements: (i) Rq = 8.33 ± 0.35 kpc 
[55] , provided by the monitoring of the orbits of S- 


stars around the central supermassive black hole, and 
(ii) = 30.26 ± 0.12km/s/kpc [43], provided 

by the proper motion of Sgr A*. Combining (i) and (ii) 
with the allowed values of Vq reported above, we have 
the following ler configurations: 

(a) Rq = 7.98 kpc, Vo = 214.52 km/s for Vq = 26 km/s; 

(b) Rg = 7.98 kpc, vq = 237.18 km/s for Vq = 5.25 km/s; 

(c) Rq = 8.68 kpc, Vo = 235.62 km/s for Vq = 26 km/s; 

(d) Rq = 8.68 kpc, Vq = 258.45 km/s for Vq = 5.25 km/s. 

We use configurations (a) through (d) to quantify the im¬ 
pact of current uncertainties on our results. This is shown 
explicitly in the bottom panels of Fig. [2] for fiducial bary¬ 
onic model I 19, 21, 31 . It is possible to conclude that, 
within the same baryonic model, the best-fit values of ao 
are relatively insensitive to the actual local circular ve¬ 
locity vo provided the Sgr A* constraint is met (cf. (a) 
and (b)). Moreover, a value of Ro on the high end of the 
currently allowed range shifts the favoured ao to higher 
accelerations, i.e. away from the range inferred from ex¬ 
ternal galaxies (cf. (c) and (d)). These considerations 
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FIG. 3: The MOND interpolating function as inferred directly from the rotation curve of our Galaxy. The red and blue 
bars represent the binned la measurement of p(a/ao) = ajv/a for fiducial baryonic models I |1T)[ ITT! I?!] and II [211 1271131 1. 
respectively, while the grey band encompasses the lcr measurements of all baryonic models. The slight shift in a of the red and 
blue bars is for visualisation purposes only. Overplotted are the standard and simple interpolating functions for cio = 10 -10 m/s 2 . 
For reference, we also indicate the values of a corresponding to Rq, R cu t = 2.5 kpc and R cu t = 4.5 kpc in the case of a flat 
rotation curve v c = vq. We have adopted here Ro = 8 kpc, Vo = 230km/s and (U, V, W) Q = (11.10,12.24, 7.25) km/s |201 . 


hold for both fi st( j and Therefore, the conclusions 

drawn in the previous paragraphs do not change qualita¬ 
tively when varying the Galatic fundamental parameters. 
The same applies when pushing the radius cut to 4.5 kpc, 
as also shown in the bottom panels of Fig. [2] 

In Fig. [3] we perform a data-driven test for MOND- 
like theories that can be applied directly to rotation 
curve and photometric data. By mapping the modi¬ 
fied acceleration with the rotation curve data through 
a = RiOc and the Newtonian acceleration with the bary¬ 
onic distribution through un = Rw 2 , we reconstruct /r 
directly from Eq. 0 with no dependence on ao, namely 
/i(a/ao) = a^r/a. Then, the resulting measurements of 
/r are binned in linear intervals of size A a = 10 -10 m/s 
requiring at least five rotation curve measurements per 
bin. This is shown in Fig. [3] by the red bars for fidu¬ 
cial baryonic model i pa nn eu, by the blue bars for 
an additional baryonic model II [2T| [27} [3Tj and by the 
grey band for the bracketing of all baryonic models imple¬ 
mented in our analysis. In the same plot we superimpose 
Mstd(a/ao) and fi s i m (a/ao) for do = 10~ lo m/s 2 . As clear 
from the figure, the standard interpolating function is in 
tension with the data for both fiducial baryonic models, 
whereas the simple interpolating function remains viable 
for fiducial baryonic model I. Note that for a flat rota¬ 
tion curve v c = vq the acceleration reads a = Vq/R, so 
the kinematic tracers at small R correspond to large a 


and vice-versa. Therefore, the importance of baryonic 
modelling grows for large a since that corresponds to the 
region of the Galaxy where the baryons become more and 
more important. Although at present uncertainties are 
still sizeable, Fig. [3] provides a powerful non-parametric 
test of the MOND paradigm. Our analysis is complemen¬ 
tary to that of Ref. BE which covers a slightly different 
acceleration range, and that of Ref. m, which is more 
focussed on the outer Galaxy. 

From Figs. [2] and [3j it appears that the simple interpo¬ 
lating function provides a better fit to Milky Way data, 
in agreement with earlier studies which adopt a different 
setup and approach, e.g. mm- Let us recall that this 
functional form is strongly disfavoured by solar system 
tests and it has to be modified in the large acceleration 
limit (see e.g. mi eh])- Modifications of the simple in¬ 
terpolating function have been proposed in the literature 
(i.e., the family of I'-functions mi or the “improved sim¬ 
ple” function m) by making /i converge to unity for 
a/a 0 > 10- We note that, for the range of accelerations 
to which we are sensitive (cf. Fig. [ 3 ]), these modified func¬ 
tions and the simple function are virtually indistinguish¬ 
able for typical values of ao and the results presented here 
for fjysim would also apply to those modified functions. As 
our measurements of the Galactic rotation curve and our 
understanding of the distribution of baryons improve, we 
expect these tests to provide more and more stringent 
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constraints on both fj, and do- For the time being, how¬ 
ever, we can conservatively state - based on an up-to-date 
compilation of rotation curve data and a comprehensive 
collection of data-inferred morphologies of the baryonic 
component - that MOND variants employing the stan¬ 
dard interpolating function do not fit simultaneously the 
rotation curves of external galaxies and that of the Milky 
Way. 

Let us finally point out that we have adopted in our 
analysis a wide compilation of baryonic models to bracket 
the current uncertainty on the morphology and compo¬ 
sition of our Galaxy. In the coming years, astronomical 
surveys such as Gaia [35], APOGEE-2 (SDSS-IV) [50] . 
WFIRST [ST], WEAVE [52] and 4MOST [S3] will even¬ 
tually help narrow down on the baryonic distribution, 
thus providing a decisive step towards testing the MOND 
paradigm in the Milky Way. 
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